Minimal field requirement in precessional magnetization switching 
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We investigate the minimal field strength in precessional magnetization switching using the 
Landau-Lifshitz-Gilbert equation in under-critically damped systems. It is shown that precessional 
switching occurs when localized trajectories in phase space become unlocalized upon application of 
field pulses. By studying the evolution of the phase space, we obtain the analytical expression of the 
critical switching field in the limit of small damping for a magnetic object with biaxial anisotropy. 
We also calculate the switching times for the zero damping situation. We show that applying field 
along the medium axis is good for both small field and fast switching times. 

PACS numbers: 75.60. Jk,75. 75. +a 



I. INTRODUCTION 

Magnetization reversal in magnetic particles and thin 
films has been a continuous and growing topic in the 
past several decades, motivated by its great application 
potential in magnetic data storage and random access 
memories (RAM). Recent development of new fabrication 
techniques has made it possible to produce nanometer- 
sized magnetic objects with well-controlled shape, struc- 
ture, and chemical composition. Magnetic anisotropy of 
these objects, including shape anisotropy and magneto- 
crystalline anisotropy, makes the dynamic magnetization 
processes highly nonlinear. A thorough understanding of 
micromagnetic dynamics is thus desirable as many efforts 
have already been made in this field. ^'^ 

In the past, magnetization reversal is realized by ap- 
plying a magnetic field pulse mainly antiparallel to the 
initial magnetization. The magnetization will then un- 
dergo multiple rotations around the local effective field 
to reach the final equilibrium direction. This process 
is best viewed by drawing the energy landscape of the 
system. Assume magnetization in the x and —x direc- 
tions gives the local energy minima. As shown in Fig. ^ 
upon application of the field pulse, the landscape is tilted, 
making the barrier disappear. The system will jump to 
the global minimum and the magnetization is therefore 
switched from one direction to another. It is clear that 
the applied field has to be strong enough to overcome the 
energy barrier. The critical switching field is first dis- 
cussed in the uniaxial anisotropy case in the pioneering 
work by Stoner and Wahlfarth,'^ and recently extended 
to non-uniaxial and three dimensional cases.^il Energy 
dissipation is necessary in this process for the system to 
move from one minimum to another. We call this type of 
reversal the Stoner- Wahlfarth (SW) type. Typical rever- 
sal time for such a process is of the order of nanoseconds. 

Recently, a novel approach towards ultrafast magne- 
tization switching by precessional motion has been pro- 
posed and observed experimentallyi.SiL2i2ii24iiii^ Instead 
of applying a field pulse antiparallel to the initial magne- 
tization, a perpendicular field pulse is applied and in- 
duces a large angle precession. If the pulse is termi- 



nated at 180° angle of the precession, whose period is 
usually of the order of picoseconds, the magnetization 
is reversed. However, the effective reversal times could 
be several nanoseconds due to the decay time of resid- 
ual magnetic precession ( "ringing" ) if the magnetization 
does not end up in the final equilibrium state precisely. 
One can overcome this difficulty by fine tuning the pulse 
parameters so that the magnetization will move along 
the so-called ballistic trajectory, eliminating the ringing 
effect. ^'^« Thereby the fundamental ultrafast limit of field 
induced magnetization reversal is reached. In contrast 
to the SW type reversal, the precessional switching is so 
fast that dissipative effects can be neglected. In other 
words, the system energy is conserved and the problem 
can be formulated using the Hamiltonian equations. 

Such a precessional switching opens a way to re- 
duce not only the reversal time but also the field re- 
quired for switching. As pointed out by early numer- 
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FIG. 1: Typical energy landscapes: energy E vs magnetiza- 
tion M at different value of applied field H. In (a) H — 0, 
system sits on one of its local minimum 1 separated with 2 by 
a barrier. In (b) H < He, the critical field, position 1 is still 
a local minimum and magnetization reversal will not occur. 
In (c) H — He, the minimum 1 and the barrier coalesce, any 
fluctuation will cause the system jump to minimum 2. In (d) 
H > He, there is only one minimum, the magnetization is 
reversed. 
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ical calculationSfi^iii precessional magnetization switch- 
ing can be observed well below the Stoner-Wohlfart limit. 
This result is recently recovered by considering the pre- 
cessional switching as a result of bifurcation, the long- 
term behavior of the dissipative system.^ However, to 
our knowledge, no systematic study of the minimal field 
strength exists in the literature and most of the work was 
done by numerical means. 

In this paper we investigate the minimal field strength 
required for precessional switching of a magnetic object 
with biaxial anisotropy. We use phase space analysis as 
our main tool to tackle this problem. It turns out that 
the mechanism of the precessional switching is directly 
related to the evolution of the localized and unlocalized 
trajectories in the phase space. Unlike the SW type 
switching, where the fixed points and their movement are 
the central objects to study, in the precessional switching 
the states move along equienergy curves and we need fol- 
low their motion globally. The minimal (critical) field is 
obtained when the localized trajectories become unlocal- 
ized. By studying the phase space evolution we are able 
to obtain the analytical expression for the critical field. 
This result provides useful information on the design of 
field pulses for precessional switching based devices. 

The paper is organized as follows. In Sec.^we present 
a simple model describing magnetization dynamics and 
define the central problem. We then in Sec. IIIII consider 
the zero damping situation in which the switching field is 
applied perpendicular to the easy axis. We also calculate 
the switching times for this configuration. In Sec. IIVI 
with the assumption of very small damping we study the 
applied field perpendicular to the hard axis configuration 
and compare our result with the standard SW model. 
Finally the paper is summarized in Sec. 



II. MODEL 

The simplest micromagnetic model is the so-called 
macro-spin model, i.e., the magnetization is uniform and 
displays collective dynamics. This model is valid for 
small object size (reaching the single domain limit) and 
low magnetic coercivity. The motion of magnetization 
M with a phenomenological damping is governed by the 
Landau-Lifshitz-Gilbert (LLG) equation. This equation 
will be used in the following dimensionlcss form: 



dm, dm 
— — ~ — m X n-cff + am x — — 
at at 



(1) 



Here m = M/M^ is the magnetization unit vector, 
hcs — Hcff/Ms is the scaled effective field, time is mea- 
sured in units of {jMs)~^, Mg is the saturation magneti- 
zation, 7 is the absolute value of the gyromagnetic ratio, 
and a is the dimensionless damping parameter. 

We consider a very thin film in the x-y plane with 
an in-plane uniaxial anisotropy, taking x the easy axis 
of magnetization. The demagnetizing field factors are 



practically equal to zero in the film plane and one per- 
pendicular to the film plane, respectively. The magnetic 
energy density is written 



w(m, h) = -ml — -Kml — h 



m 



(2) 



where K > accounts for the scaled in-plane a;-axis 
anisotropy )and h denotes applied field. Usually the 
value of K is about 0.01 for thin films because of the 
large demagnetizing field in the z direction. However, af- 
ter proper scaling the above expression can be also used 
to describe the energy density of small magnetic particles 
with biaxial anisotropy as well. In following discussions 
we assume K > 0. 

Let us consider the energy dissipation rate. Recall that 
hcff = —dw/dm, we find out the dissipation rate has the 
following form 



dw 
dF 



dw 
dm 



dm 
"dT 



1 



Ih-off X ml 



(3) 



For small damping parameter a, short pulse period, or 
small applied field, the damping term in the LLG equa- 
tion can be neglected. We note that in real experi- 
ments the damping parameter is usually of the order of 
10~^, and the field pulse duration can be adjusted to 
hundreds of picoseconds, thus providing a very good ex- 
amples of this zero damping model. The dynamics is 
described by the following equation: 



dm 

— - = -m X hcs ■ 
ot 



(4) 



In terms of the polar angle 9 and the azimuthal angle 
the energy density reads 



1 



1 



K SIT? 6 cos^ I 



2 2 " (5) 

hx sin 9 cos (j) — hy sin 9 sin (j) ~ hz cos 9 , 



and the equations of motion are 
• 1 dw 



smf 



1 dw 
sin 9 d9 



(6) 



For simplicity we consider field pulses with zero rise 
and fall times. That is, the applied field remains constant 
for a pulse duration of T. We then define the problem 
as follows. Assume initially the magnetization is along 
the -l-i direction a.t 9 = tt/2, cf) — 0, denoted by mo. 
By applying a field pulse we expect the magnetization to 
switch to the —x direction at 9 = 7r/2, = tt, denoted 
by mi. In order to achieve this both field strength and 
pulse duration have to satisfy certain requirement. Here 
we will mainly consider the field strength requirement. 



III. ZERO DAMPING: SWITCHING FIELD IN 
y-z PLANE 

We first study the a = case. In this case, the magne- 
tization cannot relax to a local minimum by dissipating 
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FIG. 2: The phase space at if — 0.2 for different /i^. Point 
X denotes the initial state mo- In (a) = 0, mo is a local 
minimum. For < hz < 1, mo gradually moves out from the 
localized region as shown in (b)-(d). In (e) 1 < /iz < 1 + K, 
the two saddle points collide at the north pole, end up with 
only one saddle point. In (f) all trajectories are unlocalized 
in the range hz > 1 + K. 



energy, therefore the initial state mo and the final state 
mi must be connected by a constant energy curve during 
the application of the field pulse. This requirement leads 
to h ■ (mi — mo) = 0; the applied field must be in the 
y-z plane. 

The phase space without applied field is shown in 
Fig. [2Jl. The system has (a) two minima at — 7r/2, 
(j) = 0, TT, (b) two maxima at 6 ^ 0,tt {(f> has no defini- 
tion), and (c) two saddle points at 6* = 7r/2, (p = ±7r/2. 
The phase space is divided by the separatrix (trajecto- 
ries that join the saddle points) into four regions. Inside 
the separatrix there are two regions, where the motion is 
periodic and localized around the corresponding minima. 
In a dissipative system these two regions are called basin 
of attraction for these two minima. Outside the separa- 
trix, in the top and bottom, the motion is unlocalized, 
i.e., the magnetization can move from the +x direction to 
the —X direction. Generally speaking, in a precessional 
switching the function of the applied field is to tilt the 
phase space structure and drag the initial state mo out 
of the localized region so the magnetization can precess 
globally. 

To better illustrate this idea, we consider two special 
cases in which the applied field is along the z direction 
and y direction, respectively. 




FIG. 3: The phase space ai K — 5 for different h^. Point X 
denotes initial state mo. In (a) 1 < hz < {l + K)/2, mo stays 
in the localized region. In (b) (f + K)/2 < hz < 1 + K , mo 
is on a unlocalized trajectory. 



In case (i), the field h is applied along the z direction. 
The phase space at X < 1 is shown in Fig. |2b-f. In the 
range oi < 1, the phase space has the same topolog- 
ical structure as = 0. There are (a) two minima at 
9 = arccos[/iz/(l + K)], (j) = 0,7r, (b) two maxima at 
6 = 0,7r, and (c) two saddle points at 9 = arccos/iz, 
(f) = ±7r/2. In the range < hz < \/K, the initial state 
mp stays in the central localized region, no magnetiza- 
tion switching will occur (Fig. I^b). If — VK, mo is 
on the separatrix and it could evolve into the other local- 
ized region, that is when the switching starts (Fig. O;). 
Note that on the separatrix the motion is however very 
slow when approaching the saddle points.. As shown in 
Fig. 123, in the range 1 < < 1 + K, the two sad- 
dle points at 6 = arccosft-z, (f) — ±7r/2 disappear, and 
9 = becomes a saddle point, i.e., the separatrix passes 
the north pole. Finally for /i^ > 1 + K there are only 
one minimum at 9 = and one maximum at 9 = n, all 
trajectories are unlocalized. 

li K > 1, i.e., the anisotropy of the easy axis is bigger 
than that of the hard axis, the two wells around the en- 
ergy minima are so deep that even after hz passes 1 the 
initial state mo still stays in the localized region, shown 
in Fig. calculation shows the minimal field strength 
that causes switching is hz = {l + K)/2, bigger than ^/K. 

Now let us consider the case where the field h is applied 
along the haty direction. The hy field breaks the y axis 
symmetry. As a consequence, the two saddle points have 
different energies. The phase space at /sT < 1 is shown in 
Fig. El In the range oiO < hy < K, there are (a) two min- 
ima at = 7r/2, </) = arcsin /ly/i^, TT — arcsin hy/K, (b) 



two maxima at 9 
and (c) two saddle points at 9 



arcsin hy,iT 



arcsin h 



J/' 



±7r/2. As 



hy increases, the size of localized regions shrinks. The 
precessional switching occurs at hy = K/2 (Fig. ^p). 
Fovhy > K, the two minima collide together with one of 
the saddle points and disappear, a new minima emerges 
at 6* = 7r/2, (p = tt/2. After hy passes 1 there are only 
one minimum at 9 = 7r/2, </> = 7r/2 and one maximum 
at 9 = Tr/2, ip = — 7r/2. If > 1, evolution of the 
phase space will be different, but the critical condition 
hy = K/2 stiU holds. 
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FIG. 5: Critical switching field in the y-z plane at a = for 
{a) K < 1 and (b) K > 1. The intersections of the curve are 
in (a) hy = K/2, h, = y/K, in (b) hy = K/2, h, = {K + l)/2. 
Fields shown in dashed line corresponds to the situation where 
the initial state has the same energy with one of the saddle 
points but is not on the separatrix. 



FIG. 4: The phase space a.t K = 0.2 for different hy. Point 
X denotes the initial state mo. In (a) and (b) < hy < K . 
At hy = K/2, mo is on the separatrix. In (c) K < hy < 1. 
In (d) hy > 1. 



In following discussion we assume this condition al- 
ways holds since we are interested in the minimal field 
strength. We will see later the result is consistent with 
this assumption. The energy of the fixed points is 



From the above discussion, we can see that preces- 
sional switching occurs during the transition of trajec- 
tories around the initial state mo from localized to un- 
localized. In other words, this is when mo moves from 
inside to outside the separatrix, so the critical field is the 
one that makes mo on the separatrix. Since points in 
the phase space move along trajectories of constant en- 
ergy, the critical condition can be also stated as that the 
energy of initial state equals the energy of saddle points 
(but we need check if they are on the same trajectory). 
We observe that this situation usually happens before 
collision of two fixed points, which is the critical condi- 
tion for the SW type reversal. This explains the smaller 
switching field in precessional switching. 

Now we consider the general case in which the field h 
is applied in the y-z plane at an arbitrary angle. The 
energies of the initial state mo and final state mi are 
always —K/2. We are looking for saddle points with the 
same energy. The fixed points are given by 



sin 9 cos 6(1 + K cos 



hy cos 9 sin < 



K sin^ 6 sin (f> cos (f> — hy sin 6 cos ( 



hz sin 9 = 0, 
(7a) 

= . (7b) 



From Eq. I|7b|l we see that solutions can be grouped 
into two categories. The first category includes solutions 
with s'm4> = hy/{K sin 9). Inserting this expression into 
Eq. H7a|l and solving it gives us fixed points 



cosf = 



1 + K ' 
with the condition 



smc 



hi 



K^l-hl/{l + Kf 



< 1 



(8) 



(9) 



E 



K 

y 



2\K 1 + K 



(10) 



The energy is always smaller than —K/2 so they are not 
what we are looking for. Stability analysis shows in fact 
they are energy minima. 

Now we turn to the other category, which includes so- 
lutions with COS0 = 0. In this category, fixed points 
occur on a great circle parametrized by setting = 
and letting 9 run from to 2tt. This makes the problem 
essentially a two-dimensional problem in the y-z plane. 

Instead of numerics, this problem can be solved by 
a geometrical approach, first proposed by Slonczewski 
in the Stoner-Wohlfarth problem, then extended by 
Thiaville^ii^ In essence, it considers the field h, instead of 
the magnetization direction m as the main variable. We 
write m = s'm9y -\- cos 9z = {sin 9, cos 9) and its orthog- 
onal vector e = {cos 9, — sin9). The energy requirement 
is 



E 



1 



■ cos 9 — h ■ m 



K 
~2 



The extremum condition is 



d^ 



sin 9 cos 9 — h ■ e = 



(11) 



(12) 



Since the system already has two minima at points ||SJ|, 
the saddle point we are looking for has to be the minimum 
along the 9 direction. This leads to 



d^E 
d6i2 



sm 



h m > 



(13) 



The critical switching field is determined by the inter- 
sections of the two lines defined in Eqs. (|ll|l and l|12|l . 
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FIG. 6: Switching times contour in the y-z plane for K = 0.2. 
The dashed line in the center is the critical switching field. 
From outside to inside, the switching times are 2, 3, 4, 7, 
10, 20, measured in {'yMs)~^ . For typical configuration it is 
about 0.01 ns. 



We find 

hy = ^ smO{K — cos^ 6) , 

hz = ^ cos ^(l + A' + sin^ 9) . 
Inserting above expressions into Eq. (|13|) yields 



d6l2 



= 1 



-K 



> . 



(14a) 
(14b) 

(15) 



For K > 1, the abov e inequality always holds. For K < 1 
we have |cos6'| < + K) /3. However, in this range 
the critical curve does not form a closed loop, shown in 
Fig. The reason is that at nonzero hy the two saddle 
points have different energies so the trajectory starting 
from mo may not pass that saddle point. After check- 
ing we can remove these "faked" solutions by requiring 
Icos^l < Vk. 

From Fig. O we can see that if if ^ 1 the magnitude 
of critical switching field along any direction in the y-z 
plane is about the same; for magnetic thin films where 
K <^ 1 then y direction is the best choice to apply small 
field. We also calculate the switching time for K < 1 hy 
numerically integrating Eq. ((HJ. The equi-timc contours 
are plotted in Fig.El As the magnitude of fields increases, 
the switching time decreases quickly from the center. 



IV. SMALL DAMPING: SWITCHING FIELD IN 
x-y PLANE 

In this section we study the small damping case. As 
finding the exact analytical solution to the equation of 



magnetization motion including damping is a formidable 
task, we turn to study the limiting case, where a is small 
enough so that we can still treat the problem as a Hamil- 
tonian problem and a is also big enough so that after 
the application of the field pulse the magnetization can 
relax to the local minimum in finite time. Because of 
damping the energies of the initial and final state need 
not to be the same, a nonzero field component along the 
X direction is allowed. 

We choose to apply the field in the x-y plane. By this 
way, we can compare our result with the SW model. The 
basic procedure is the same as last section: we will search 
the saddle points with the same energy of the initial state, 
then map to the whole switching field. 

It can be seen from Fig.^that the two fixed points that 
originate from the south and north pole do not involve 
in this precessional switching, so the saddle point we are 
looking for is on the x-y plane too. Let 9 = 7r/2, then the 
problem turns into a two-dimensional problem in the x-y 
plane. We write m — cos (jjx H- sin (jjy — (cos 0, sin </)) and 
its orthogonal vector n = (— sin 0, cos (/>). The applied 
field is h = {hx,hy). 

The energy requirement is 



E 



K 



' — h ■ m — 



K 



The extremum condition is 
AE 



—— — K cos (h sin 6 — h ■ n = Q 



(16) 



(17) 



Solving above equations gives us the switching field in 
the x-y plane: 



hx 



K 



C0S( 



K 



■ sin Ml 



(18a) 
(18b) 



In order to check the stability of the corresponding fixed 
point we carry out the standard routine. With h^, hy 
taking the above form, the second derivatives of w are 



d9d(p 



= 1 



K 



(1 — cos < 



(19a) 



A:(cos0- l)(cos(^+ -) , (19b) 



= 



(19c) 



It is obvious that has to in the range < (/) < |7r 
(upper branch in Fig. [TJ or— |7r<0<O (lower branch 
in Fig.dl to make the fixed point a saddle point. Again 
we still need to check if the saddle point and the initial 
state are on the same trajectory. We require both hx and 
hy are monotonous function of in the upper and lower 
branches, which narrows the range to — f'"' < 4> < |7r. 

To compare with the switching field in the SW model, 
which takes the form h = {— cos^ (j), sin^ (f>), we plot both 
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FIG. 7: Critical switching field in the x-y plane for switching 
from the x direction to the —x direction. The solid line shows 
the switching field of precessional switching at zero damping 
limit. The dashed line is the switching field of the Stoner- 
Wohlfarth model. For completion we also draw the right half 
of the SW astroid. Arrow indicates the minimal field strength. 

fields in Fig.d It is clear that precessional switching can 
occur well below the SW limit. Surprisingly, magnetic 
switching is possible for a field with a positive x compo- 



nent. Although in this case the minimal field strength 
is 0.38K (coscj) = 1/3), considering the long relaxation 
time it is still better to apply the field mainly along y 
direction with a small x component to compensate the 
energy loss during the motion. 



V. SUMMARY 

The critical switching field of a monodomain magnetic 
object with biaxial anisotropy has been studied. In par- 
ticular, switching fields perpendicular to the easy and 
hard axes are calculated using the geometrical method. 
It is found that in the limit of zero damping, applying 
field along the medium axis is good for both small fields 
and fast switching times. 

This method can be generalized for magnetic parti- 
cles with a higher order anisotropy if one follows the 
same steps outlined in this paper. The two-dimensional 
model may not be sufficient because one needs to consider 
more fixed points. A more complicated three-dimensional 
model would then be required. 
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